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Abstract 

Bright bursts of gamma rays from outer space have been puzzling Astronomers 
for more than thirty years and there is still no conceptually complete model for the 
phenomenon within the standard model of particle physics. Is it time to consider a 
supersymmetric (SUSY) origin for these bursts to add to the astronomical indications 
of supersymmetry from dark matter? 



1 Introduction 



Several times each day, satellites in near earth orbit observe a prodigious burst of gamma 
rays and the accumulating data has resisted a complete physical explanation for decades. 
How long should one wait until deciding that it is time to consider explanations beyond the 
standard model? Astronomers, while recognizing a persistent mystery behind the "central 
engine" of violent astrophysical events, have been notoriously slow to encourage discussion 
of new physics solutions. 

Ideally, in a situation like that of the gamma ray bursts (grb's), there should be a broad 
search for conceptually complete explanations that predict the primary characteristics of the 
phenomenon in zeroth order. Violent astrophysical events are complicated processes and one 
must expect difficult numerical computations to be necessary even if the underlying physical 
model is conceptually complete. By conceptually complete I mean that each step of the 
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process is completely based on known or plausibly proposed physical principles. We suggest 
that supersymmetry does provide a conceptually complete picture of the burst phenomenon 
while the standard model of particle and astrophysics does not. 

Within the standard approaches to gamma ray bursts there are significant conceptual 
gaps. Among these are the mechanisms for the enormous energy release required, the mech- 
anism for the transfer of this enormous energy into a narrow range of the gamma ray spec- 
trum, and the mechanism for the strong angular collimation of the bursts if such exists or 
the extra energy release otherwise required. The places where more physics input is needed 
are currently marked by the signs "central engine" , "firecones" , "sub-jets" etc. 

In addition, the standard approaches do not predict the primary quantitative character- 
istics of the gamma ray bursts except as related to free parameters in the theory. Most of 
the studies of these bursts take as a starting point the unexplained production of relativistic 
outgoing jets by a dying star and proceed to model the afterglow left in their wake. Several 
good reviews PQ of the standard approaches are available. 

The primary characteristics of the phenomenon are: 

1. Dominant photon energies in the E~ — 0.1 to 1 MeV range. 

2. The distribution of burst durations ranges from approximately 0.02 s to 300 s with 
a pronounced dip [3] in the distribution at 2 s. Bursts of less than 2 s duration are 
referred to as "short" bursts and those above 2 s are defined to be "long" bursts. 

3. Total burst energy of the order 



where AQ is the opening solid angle of the burst. If AQ is Att (isotropic burst), 
the energy output is equivalent to or greater than that observed in supernovae. The 
difference, of course, is that supernovae put most of their energy release into neutrinos 
and into kinetic energy of heavy particles and nuclei, eventually appearing in a broad 
range of the electromagnetic spectrum. Subject to some possibly strong assumptions, 
arguments have been given jl] that the jet opening angles are ~ 5 deg. This would 
correspond to AQ « 7.6 • 10~ 3 . and E? ~ 1.8 • 10 50 ergs. For comparison with the 
implied burst energies, the rest mass of the sun is 1.9 • 10 54 ergs and that of the earth 
is 6 ■ 10 48 ergs. 

Thus, the shortest bursts can be compared to a hypothetical near solar mass of matter 
disappearing into MeV scale photons in 20 ms. Relativity would then require that this solar 
mass should be initially compressed into the radius of the earth (Re/c = .02s) or less. Such 
is the typical density of white dwarf stars. However, the sun is prevented from converting 
its full rest energy to photons by baryon number conservation which is well tested on time 
scales up to 10 32 yr. In addition, compact objects such as white dwarfs are Fermi degenerate 
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systems restricted in the energy they can emit by the Pauli principle. The electrons occupy 
all energy levels up to the Fermi energy. The average and total kinetic energy of a degenerate 
electron gas in a white dwarf of solar mass and earth radius are, respectively 0.1 MeV and 
1.2 • 10 50 ergs, not far from the observations quoted above assuming a 5° opening angle. 

This suggests that one consider the elimination of the Pauli blocking by a transition to 
a system of Bosons as is possible (perhaps even necessary) given supersymmetry and other 
physical considerations A neutral system of density p, fully relieved of Pauli blocking, 
will collapse to a black hole in a time 



Thus one might expect that the bursts of duration somewhat above 1.5 s will be depleted 
and might cause a dip in the duration distribution as is observed at 2 s. Bursts of duration 
much less or much longer than this time might be relatively unaffected by the cutoff due 
to gravitational collapse. A full treatment, however, requires studying the collapse during a 
gradual lifting of the degeneracy pressure. 

Although a SUSY phase transition to a system of Bosons can easily explain the primary 
characteristics of the gamma ray bursts as outlined below, my astronomer friends tell me 
that substantial progress is also being made in the standard astrophysical approaches to the 
short duration bursts. Among the most recent ideas that have been proposed is the following 



Depending on various parameters, up to 10 52 ergs of energy could be released into a 
neutrino-antineutrino cloud which could then be converted on a short time scale into an 
e + e~ cloud and made available for the production of a relativistic fireball. Assuming a one 
percent conversion efficiency, the fireball would have about the right rest energy, 10 50 ergs, 
to account for the gamma ray bursts. It would seem that such standard model explanations 
for the bursts still do not have the predictive power or the conceptual completeness of the 
SUSY phase transition model. Similarly, standard model attempts at a unified model of the 
bursts with a dip in the duration distribution at 2 s often assume the dip is due to a viewing 
angle dependence [DEI- This, however, would seem to require a somewhat strange and ad 
hoc shape of the gamma ray wave packet. 

The SUSY phase transition model for gamma ray bursts is based on the following scenario. 

1. In a region of space with a high level of fermion degeneracy there is a phase transition to 
an exactly supersymmetric ground state. In the SUSY phase, electrons and their SUSY 
partners (selectrons) are degenerate in mass as are the nucleons and snucleons, photons 
and photinos etc. We assume that these masses are at most those of the standard 
model particles in our normal world of broken supersymmetry. This assumption is, 
perhaps, supported by the superstring prediction of low mass (massless) ground state 
supermultiplets. 

2. In the SUSY phase, electron pairs undergo quasi-elastic scattering to selectron pairs 




(1.2) 



e(pi) + e(p 2 ) -> e(p 3 ) + e(p 4 ) 



(1.3) 



which, uninhibited by the Pauli principle, can fall into the lowest energy state via 
gamma emission. These gamma rays are radiated into the outside (non-SUSY) world. 
Selectrons and other SUSY particles are confined to the SUSY bubble being too low 
in energy to penetrate into the broken SUSY world where their mass would be orders 
of magnitude greater. Since there are about 

iV = 6-10 56 (1.4) 

electrons in a white dwarf and the average squared momentum in the degenerate sea 
is 

Q / Q AT \ 2 / 3 

<P 2 >=1 (^-y) ■ (2^) 2 « 0.15(MeU/c) 2 , (1.5) 
the total energy released if all the electrons convert is 

E T = N(< E > -mc 2 ) = 1.2 • 10 50 MeV (1.6) 

as advertised above. 

3. The highly collimated jet structure might be produced by the stimulated emission of 
sfermions and photons; in a bath of pre-existing sfermions, new sfermion production 
goes preferentially into momentum space bins with large occupation numbers. 

4. Simultaneous with electron conversion into selectrons, nucleons within heavy nuclei 
convert into snucleons. With no further support from the electron degeneracy, the star 
collapses to nuclear density under gravitational pressure. 

5. Remaining nucleon pairs then undergo the analogous conversion to snucleon pairs with 
the cross section mediated by the strong exchange of supersymmetric pions. This 
process can be temporarily interrupted by brief periods of fusion energy release but 
then continues until all kinetic energy is radiated away or until the star falls below the 
Schwarzschild radius and becomes a black hole. 



In this model many of the bursts are due to the decay of isolated white dwarfs or neutron 
stars which are absolutely table in standard astrophysics. We therefore predict the existence 
of low mass black holes below the Chandrasekhar limit. Other bursts might come from the 
SUSY transition during a stellar contraction when a stage of fermion degeneracy is reached. 
If the duration is linearly related to the radius of the star, the longest duration bursts would 
come from a star of radius 10 4 times that of earth (10 2 times that of the sun). However, if 
the photons undergo a random walk with a mean free path of earth radius before emerging 
from the dense star, durations would depend quadratically on radius 
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In such a case, the longest duration bursts would come from stars of radius 100 times that of 
earth. Many other effects, however, including fusion reignition and the free collapse time of 




Figure 1: Effective potential in vacuum. The broken SUSY world is indicated with a postive 
vacuum energy density while the SUSY phase has zero ground state energy density. 

a star relieved of Pauli blocking, can also influence the duration of the burst so the situation 
needs more study. 

There are now accumulating indications from accelerators and astrophysics that we live 
in a broken SUSY world with a positive vacuum energy e « 3560 MeV/m 3 . On the other 
hand string theory seems to persistently predict that the true vacuum is an exact SUSY 
world with zero vacuum energy. In this situation it is inevitable that bubbles of true vacuum 
will spontaneously form in our world and, if of greater than some critical size, R c , will expand 
to engulf the universe. The potential of some effective scalar field may be as indicated in 
figure ^ where our local minimum of broken SUSY is labelled "US" and the true vacuum 
is exactly super symmetric. In reality, $ may stand for a complicated linear combination of 
string moduli and there may be a large number of local minima most of which are of a height 
(D(10 100 ) times greater than e. 

In the simplified single field situation the critical radius is obtained by an instanton 
calculation [Sj. 

R c =-. (1.8) 

e 

Here S is the surface tension of the bubble of true vacuum. Since it is to be expected 
that the probability to nucleate a bubble of radius r is a steeply falling function of r, for 
sufficiently large - the broken SUSY phase is metastable [5| in dilute matter. Many authors 
have discussed a catalysis of the phase transition in the presence of dense matter [TUj, [TTj . 
If the effective field, $, is a gravitational modulus, it is reasonable to assume that the critical 



radius is related to the difference in ground state energy densities as suggested by figure 121 
that is 

R c = (1.9) 
e + Ap 




Figure 2: Effective potential in the presence of matter. 

where Ap is the ground state matter density, p in the broken SUSY phase minus the 
corresponding ground state matter density ps in the SUSY phase. In a degenerate star Ap 
is the kinetic energy density of the Fermi gas since, in the SUSY phase, the ground state 
will consist mainly of Bosons in the lowest energy level. The phase transition to the exact 
SUSY vacuum is then enormously exothermic. For the typical white dwarf 

Ap = 6 • 10 u MeV/m 3 (1.10) 

Thus we might reasonably expect that the ratio of the critical radius in a white dwarf star 
to that in vacuum is 
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Here we have made the simplest assumption, namely that the surface tension is independent 
of the density. Other possibilities might be explored. If the critical radius [3] in a white 
dwarf is of order 1CU 4 m, then the critical radius in vacuum is many times larger than the 
galactic radius. Thus a supercritical SUSY bubble in a white dwarf would grow only to the 
edge of the star where it would then be stabilized from further growth. 



A transition to a Bosonic final state as proposed in this model not only allows the release 
of the fermion kinetic energy but also suggests significant jet structure due to the Bose 
enhancement effect familiar from laser physics. 

In addition to the momentum dependence of the elementary process, the matrix element 
for the emission of a selectron pair with momenta p 3 and p 4 in process 11.31 in the presence 
of a bath of previously emitted pairs is proportional to 

M ~< n(p 3 ) + l,n(p 4 ) + 1 | a^(p 3 )a ] (p 4 ) | n(p 3 ),n(p 4 ) > 

~ y/(n(p 3 ) + 1)VW 4 ) + !)• (1-12) 
The cross section is, therefore, proportional to (n(p 3 ) + l)(n(p 4 ) + 1) . 

The calculation [T2j of the process 11.31 has been recently extended [IBj to include the 
electron mass effects which are important in the exact SUSY phase. 




Figure 3: Feynman graphs for the conversion of an electron pair to selectrons via (Majorana) 
photino exchange. 



There are two types of scalar electron (partners of the left and right handed electron) 
which must be treated separately. The three different final states are found to have the 
elementary squared matrix elements 

\M LR \ 2 = £ ss , (M a M\ + M h M\ + M a M\ + M b M\ 
1 
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where S ss / denotes averaging over the spins of incoming electrons. 
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For example the conversion rate per unit volume in the degenerate electron gas into left 
selectrons is given by 



77- = ■ ■■ =dt?^^—dp 3 dcos(6 3 )d<j) 3 dp A dcos(6 A 

V (27r)y s (s-4m 2 ) 2E 3 2E 4 

(n L (p 3 ) + l)(n L (p 4 ) + l). (1.15) 

An event generator for this process has been constructed |T3] where, in principle, one gen- 
erates one event at a time continually updating the probability distribution in accordance 
with the Bose enhancement factors. Initially, all the occupation numbers are zero but, after 
the first event, the next event is four times as likely to be into the same phase space cell 
as into any other. Because of the huge number of available states, the second transition 
is still not likely to be into the first phase space cell. However, as soon as some moderate 
fluctuation of selectrons has been produced in a single cell, the number in that state escalates 
rapidly producing a narrow jet of selectrons. As these decay down to the ground state via 
bremstrahlung, a narrow jet of photons is created around the direction of the initial selectron 
jet. After a large number of events, a jet structure emerges as the distribution locks in on 
particular cells in the final state phase space. To speed up the emergence of the jet structure, 
the occupation numbers are, in some studies, incremented by two units at each event. The 
distribution in momentum magnitude, polar angle cosine, and azimuthal angle after 6 ■ 10 5 
events is shown in table 1. The selectron momentum will essentially be totally converted to 
gamma rays as the selectrons fall into the ground state although we have not as yet treated 
this radiative process in detail. 
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0.324 


2831 


0.30 


1272163 


4.084 


29715 
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0.70 
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0.90 


282697 
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Table 1: Selectron momentum and angular distributions showing the effect of boson en- 
hancement 



Interesting features of the distribution are the jagged nature of the angular distribution 
and preference for momenta above the average of the Fermi sea. This latter effect suggests 
that earlier (more probable) emissions will be of greater energy than later ones in agreement 
with indications from burst observations [T4*j . 



There are obviously many details of the behavior of a dense supersymmetric system re- 
maining to be investigated but the SUSY phase transition picture provides an alternative 
framework for thinking about gamma ray bursts that predicts some of the basic character- 
istics and addresses some of the conceptual gaps of the standard approaches. 
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